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We study the effect of temperature (T) and baryon density (/r) dependent hadron masses on the 
thermodynamics of hadronic matter. We use linear scaling rule in terms of constituent quark masses 
for all hadrons except for light mesons. T and p dependent constituent quark masses and the light 
mesons masses are computed using 2+1 flavor Nambu-Jona-Lasinio (NJL) model. We compute 
the thermodynamical quantities of hadronic matter within excluded volume hadron resonance gas 
model (EHRG) with these T and p dependent hadron masses. We confront the thermodynamical 
quantities with the lattice quantum chromodynamics (LQCD) at p = 0 GeV. Further, we comment 
on the effect of T and p dependent hadron masses on the transport properties near the transition 
temperature (Ik). 
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I. INTRODUCTION 

The hadron resonance gas model (HRG) and its extended version, excluded volume HRG (EHRG) is the simplest 
and the most successful effective model of quantum chromodynamics (QGD) describing the hadronic phase of strongly 
interacting matter. The HRG model is essentially based on the result that the interacting hadron resonance matter 
can be approximated by that of non-interacting gas of hadrons and all the resonances [l|. Before this, there was an 
interesting result obtained in Ref. ii where the authors showed that the interacting pion gas behaves similar to the 
ideal system containing pions and p mesons. The HRG model has also been confronted with the lattice QGD (LQGD) 
and found to be in a very good agreement^. It has been extensively used in the literature for various purposes, viz., 
to obtain hadron yield in heavy ion collision experiments [^-Q, estimating transport properties of hadronic matter like 
shear viscosity [ol-ll^ and the study of fluctuations in conserved charges in heavy ion collision experiments . 

Albeit non-interacting HRG is in agreement with LQGD, it misses one important feature of hadronic interactions; 
the repulsive interactions. The existence of repulsive interactions has been confirmed from nucleon-nucleon scattering 
experiments and, in fact, used to assign finite hard-core radius to nucleons [H, [2^. Further necessity to include short 
range repulsive interaction comes from heavy ion collision experiments. It has been observed that the chemical freeze- 
out parameters obtained from fitting the particle number ratios at AGS and SPS energies lead to large values of total 
particle number densities. One can suppress these number densities by incorporating the repulsive interactions in 
HRG model through excluded volume corrections (23 - [^ . This HRG model with excluded volume correction (EHRG) 
has been found to be in good agreement with LQGD up to temperature, T^ 0.140 GeV[^ [^. 

Since HRG (as well as EHRG) is a statistical model, the essential starting point is to find the partition function 
which in this case is just the partition function of an ideal gas summed over all the hadronic states and their resonances. 
While calculating the partition function at temperature (T), it is the zero temperature (and baryon chemical potential) 
hadron masses (M/j) that enters the Boltzmann factor, Exp(-M;j/T). It is well established fact that the chiral symmetry 
is an essential feature of QGD, the spontaneous breaking of which is responsible for the large part of the quark mass 
called constituent quark mass, whence the hadrons. Further, LQGD as well as other effective model calculations 
at finite temperature shows that this symmetry is restored above so called chiral transition temperature (Tc) above 
which masses of approximate Goldstone modes drops down to current quark mass. Thus, since hadrons are made 
of quarks whose mass depends on temperature and chemical potential, it is T (and p) dependent hadron mass that 
should enter the partition function of HRG before computing any thermodynamical quantity. As we will see, taking 
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into account this effect drastically changes the thermodynamics of hadronic matter at moderately high temperature. 
Since the HRG model has been used to calculate the transport properties as well, they are also non trivially affected. 

In this paper, we use SU(3) NJL model to compute masses of constituent quarks (u,d,s) as well as low lying mesons 
(tt, K, r] and rj'). Since masses of the heavier mesons and baryons cannot be obtained in similar way, we use linear 
scaling rule in terms of constituent quarks for them. We compute all the thermodynamical quantities in EHRG at 
zero baryon chemical potential and confront them with the recent lattice QCD results^. Further, we estimate shear 
viscosity coefficient within molecular kinetic theory generalized to gas of relativistic particles. Our basic motivation 
behind this work is to emphasize the fact that the hadron masses are temperature and density dependent which should 
be used to study the thermodynamics if the HRG model is used to estimate the same especially at temperatures where 
the effects of chiral symmetry become apparent. Further, although traditional HRG model is consistent with LQCD, 
it fails to include repulsive interactions which, we feel, are important too. It is important to note that the way we 
have included T and /r dependent hadron masses is rather crude whence subjected to further improvement. 

We organize the paper as follows. In Sec. H we briefly describe the thermodynamically consistent excluded volume 
hadron resonance gas model. In Sec. HI we present the results and make comment on the shear viscosity coefficient. 
Finally we summarize and conclude in Sec. IV. 


II. HADRON RESONANCE GAS MODEL (HRG) WITH THE MEDIUM MODIFICATION OF 

HADRON MASSES 

The central quantity required to compute the thermodynamical quantities is the partition function given as 

logZ{V,T,fi) = J dm{pb{m)logZb{V,m,T,^i) +pf{m)logZf{V,m,T,^i)) (1) 

where the gas of hadrons is contained in volume V, at temperature T and baryon chemical potential p. pb and pf are 
the spectral densities of bosons and fermions respectively. Further, Z;, and Zy are the partition functions of bosons 
and fermions respectively. The hadron properties enters into the model through spectral densities Pb/f- The (free) 
hadron resonance gas model is based on the result that the thermodynamics of interacting hadron resonance gas can 
be approximated by that of non interacting gas of hadrons provided all the resonances are included in the partition 
function[l|. For such system the mass spectrum can be assumed to be sum over discrete hadronic states. 

Ma<\ 

Pb/f = ^ gaS{m-Ma) (2) 

a 

From the partition function defined by Fq. © together with the mass spectrum given by Fq. m, the thermodynamical 
quantities can be readily obtained, viz., Pressure P{T,p) = riimy_>oolii2^(V,T,/r)/V, Baryon number density ub = 
dP{T, p)/dp, entropy density s{T, p) = dP{T, p)/dT, energy density e(T, p) = Ts{T, p) — P{T, p)+ pnsiT, p), sound 
speed Cg{T,p) = dP{T, p)/de{T, p). 

In the canonical non-interacting HRG model, one can take into account short range repulsive interactions via 
Van-der-Waals(VDW) correction in the volume of the system, i.e by substitution V—>■¥—nN in the partition function 
given by Fq. ©■ Here, v = is the proper volume parameter of the hadron with the hard core radius 

With this VDW corrected HRG model, i.e excluded volume hadron resonance gas model (EHRG), one obtains the 
transcendental equation for pressure as 


P^^{T,p)=P,deal{T,p) (3) 

where Pideal is the pressure computed within canonical HRG model and fl = fi — jj,) is the effective chemical 

potential. The number density, energy density and entropy density can be readily obtained as 


^ If we consider two hadrons with the hard core radius just touching each other, they cannot come closer than the distance 2r^ . Thus, 
for each pair of hadrons the excluded volume is Vpair = |7r(2r/i)^. Whence for each hadron the excluded volume will be half of this i.e 
V = Vpair I which tums out to be 4^7rr^ . 
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In the temperature range where we are interested, the classical Boltzmann approximation is rather good approxi¬ 
mation. In this approximation, VDW prescription is merely equivalent to the factor of Exp(—uP^^/T) to the ideal 
gas pressure. 

We further extend excluded volume HRG model by taking into T and p, dependent hadron masses. To this end 
we use 2-1-1 flavor Nambu-Jona-Lasinio (NJL) model to obtain T and /i dependent masses of light mesons ( 7 r,K, 77 , 
77 ’) [ 33 , HE, [ 13 . For general reviews on NJL model, see [33 - 1^ . We cannot treat heavier hadrons and resonances at 
par wit h lig ht mesons to obtain their medium masses. So, for our purpose we use mass formulas bases on the work of 
Leupold(40l| and its generalization by Jankowski et al.jjlj. Since the hadrons are made of either two or three quarks 
to make mesons and baryons respectively, we assume that the hadrons masses scale as 


Mh = {Nq - Ns)Mq + NsMs + Kh (7) 

where Nq is the total number quarks in a given hadron and is the parameter which measures the strangeness 
content of the hadron. For the open strange hadrons. Ns is simply number of strange ( antistrange) quarks. For 
hidden strange hadrons, Ng = 2/3 for flavor singlet state while for flavor octet state, Ng = 4/3. Further, Mq and 
Mg are the constituent masses of light (u,d) and strange (s) quarks respectively. Kh is the state dependent constant 
independent of current quark mass, niq. In writing scaling rule 0 we have assumed isospin symmetry i.e = md- 
Fq. © is used for all the hadrons except for light mesons. 

For the gas of hadrons at finite temperature and baryon chemical potential, we generalize the scaling rule 0 as 

Mh{T,fd) = {Nq - Ng)Mq{T,fi) + NgMg{T,fx) + ( 8 ) 

We separate M(T = 0, /r = 0) part in above equation and absorb kh in it to get 

Mh{T,fi) = Mh{T = 0,77 = 0) + (TV, - Ng)Mq{T, ^x) + NgMg{T, ^l) (9) 

where Mq ^ is only medium dependent part of the constituent quark mass. 


III. RESULTS AND DISCUSSION 

We use the parameter set of Ref.fs^l to compute the masses of constituent quarks and light mesons using NJL 
model. To estimate the thermodynamic quantities using FHRG, we have taken all the hadrons and their resonances 
with the mass cutoff A = 2.252 GeV for baryons and 2.011 for mesons [ 4 ^. The only unknown parameter which 
fixes the excluded volume HRG is the hardcore radius Vh or the proper volume parameter v. It is customary in the 
literature to use uniform values of hardcore radius for all the hadrons [T^ H^. Baryonic hard core radius can be 
extracted from the short range repulsive interactions in nucleon-nucleon scattering processes. While it is legitimate to 
set hard core radius of all the baryons equal, detailed information regarding short range interaction between mesons 
is absent. Nevertheless, one can set same hard core radius to all mesons as that of baryons since meson charge radii 
are similar to the baryons [l^. But for our purpose we use the mass dependent hardcore radius as in Ref. [3^. In 
this scheme of parametrization, v is chosen to be proportional to the mass of each hadron; v =M/eo, eo is a constant 
which we fix to the value 0.9 GeVfm“^. We further generalize this scheme by taking into account T and 77 dependent 
hadron masses. 
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FIG. 1: (Color online)Left panel shows constituent quark masses as a function of temperature. Right panel shows light meson 
masses as a function of temperature. 





FIG. 2: (Color online)Results for thermodynamical quantities at /r = 0 GeV in EHRG-I and EHRG-II models with the mass 
dependent excluded volume parametrization. Black dashed curve corresponds to non-interacting hadron resonance gas model. 
The lattice data (red circles) is from Ref. Q. 


Fig. dKa)] shows constituent quark mass as a function of temperature at /i = 0 GeV. We note that Mu,d decreases 
with with T and drops to current quark mass around 0.2 GeV. Albeit the strange quark mass (Ais) also decreases with 
T, it does not attain its current quark mass. The temperature at which constituent quark masses attains its current 
quark mass is called transition temperature (Tc) corresponding to chiral restoration phase transition. [Fig. [TJb)] 
shows low lying meson masses as a function of temperature at /i = 0 GeV. The masses of (approximate) Goldstone 
modes show a weak dependence on the temperature as these are related to approximate SU(3) chiral symmetry 
breaking. In fact, one can show that for low temperatures dependence of mass of e.g. pion can be written as, using 
GOR relation at finite temperature. 


1 + 
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FIG. 3: (Color online)Scaled pressure as a function of temperature at two different chemical potentials in EHRG-I and EHRG-II 
with mass dependent excluded volume parametrization and eg = 0.7 GeV fm“^. 


Beyond the critical temperature the masses of pion and its chiral partner sigma become degenerate and increase 
with temperature. On the other hand, the Ua{^) is not restored within the model and pion and eta masses are non 
degenerate. Further, cusps in Fig. [Hb)] corresponds to Mott transitionjs^, [H, [s^. 

Results of the thermodynamical quantities are shown in Figs. ©. We call EHRG without T, /i dependent hadron 
masses as EHRG-I and that with T, /r dependent hadron masses as EHRG-II. We note that the thermodynamical 
quantities computed within EHRG-II start deviating from EHRG-I at O.II GeV and this deviation is more 
pronounced above 0.14 GeV. All the thermodynamical quantities are numerically larger in EHRG-II than in 
EHRG-I. This observation can be explained by simply considering Boltzmann factor Exp(-Mft,(T,/x)/T). This factor 
is a measure of probability that the specific hadronic species of mass M is thermally excited at given temperature 
whence making a contribution to the thermodynamical quantities. Since masses of all the hadrons but pions, kaons 
and eta mesons decreases with temperature, they can be thermally excited abundantly with ease. As shown in Eig. 
Hb), masses of the (approximately) Goldstone bosons do not change much around Tc, but the constituent quark 
masses do change significantly at this temperature [Fig. [Ija)]. In fact, it drops down to its current quark mass at T^. 
Since we expressed the heavy mesons and baryons masses in terms of constituent quarks [Eq. ([7])] which contribute 
significantly at higher temperatures, we see the effect of T (and /r) dependent hadron masses on the thermodynamics 
only at higher temperatures especially around transition temperature, while this effect is small at low temperatures 
where the pions and kaons are the dominating degrees of freedom. 

We note from Fig. [7] that the thermodynamical quantities computed within conventional non-interacting hadron 
resonance gas (HRG) model better fits the lattice data than EHRG-I model. Such non-interacting HRG model 
corresponds to u = 0. Thus, the lattice data seems to prefer zero excluded volume parameter, whence the point 
particle picture of hadrons. But this observation does not invalidate EHRG-I model altogether. It has been shown 
in Ref.[32| that non-interacting HRG model is problematic since the thermodynamical quantities rises very rapidly 
with temperature and ultimately shows sign of Hagedorn divergence around Tc, while in EHRG-I thermodynamical 
quantities rise less steeply than that in free HRG. Further, the better agreement of HRG model over EHRG-I with the 
LQGD may be mere coincidence since we know from the experiments that nucleons, at least, are not point particles 
but they do have finite spatial extensionfl^. There is another experimental evidence that goes in favor of EHRG-I 
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model. The analysis of the data for particle number ratios of Au+Au (AGS)^ and Pb+Pb (SPS)^ collisions suggest 
necessity to include excluded volume corrections in free HRG mo del [^. Further, including the repulsive interactions 
via. excluded volume corrections in free HRG model and with the proper choice of excluded volume parameter, it is 
observed that EHRG-I agrees with lattice data up to 0.14 GeVjS^. Eventually, by including medium dependent 
hadron masses, the agreement of EHRG-I model with LQCD is observed at higher temperatures as may be noted 
from Fig. [2ja) where the pressure (normalized by T^) computed with medium dependent hadron masses agrees with 
lattice data up to 0.2 GeV. Fig. ([3|) shows scaled pressure at finite chemical potential. We note that the pressure 
rises more rapidly even at finite /i in EHRG-H than in EHRG-I. We compare these results with LQCD simulations at 
finite chemical potentialji^. We note that at /i = 0, the estimations of EHRG-H are in better agreement with LQCD 
simulations that that of EHRG-I. Further, although the general behavior of the scaled pressure is similar we observe 
notable difference between LQCD simulations and EHRG-H estimations at /i = 0.3 GeV. But, it needs to be noted 
that these lattice simulations are estimated up to 

Fig. HKb) shows trace anomaly (interaction measure) computed within EHRG-I and EHRG-H. We note that the trace 
anomaly rises rapidly in EHRG-H as compared to EHRG-I at high temperatures. Trace anomaly in EHRG-I shows 
decreasing behavior at high temperatures. The reason behind this is twofold. First, the suppression factor ;j) 

which decreases as temperature and chemical potential increases and hence all the thermodynamical quantities in 
EHRG-I are numerically smaller than that in ideal HRG. Since there is no such suppression factor in HRG trace 
anomaly rises monotonically. Further, due to finite size of hadrons the pressure of hadron gas rises more rapidly 
as compare to energy density whence the interaction measure decreases at high temperature in EHRG-I. Strong 
suppression of thermodynamical quantities has also been observed earlier in Ref. |45l | where the authors studied 
EHRG with uniform hard core radius for all the hadrons. Although we have used different scheme of parametrization 
for hard core radius, suppression effect is still there. But in case of EHRG-H, since hard core radius is itself depend 
on temperature, the suppression effect is somehow diluted. Rapid rise of trace anomaly has also been observed in 
HRG model as well as extended HRG model which include continuum spectrum of hadrons (Hagedorn states) along 
with discrete snectrumfl^. 

Although our main purpose of this study is not to fit the lattice data, this observation is rather crucial because 
as mentioned earlier, EHRG-I fails to explain the lattice data above T^ 0.14 GeVfs^. In our previous work in Ref. 
[13, we studied the extension of HRG model by including the Hagedorn density of states at finite temperature and 
density. We found rather good agreement with the lattice data of Ref.ji^ below T= 0.15 GeV. In Ref. [i^ authors 
studied the two extensions of HRG model, viz., HRG model with excluded volume effects (EHRG) and HRG with 
continuum mass spectrum (Hagedorn states) along with discrete mass spectrum of the hadrons. They observed that 
two models are not in agreement with the lattice QCD results of Ref. when considered separately. But when 
considered together, the suppression effects in EHRG and the enhancement effects due to Hagedorn states in HRG 
leads to better agreement with LQCD. From our observation that merely including the medium effects of hadrons in 
EHRG fit the LQCD quit well, it may be tempting to conclude that the effects of Hagedorn states can be alternatively 
simulated by including T and /r dependent hadron masses in EHRG. 

The shear viscosity of the gas of relativistic particles of hard sphere radius can be calculated as [13 


T] = 


64\/8r 


■E' 


Pa 


( 11 ) 


where (| Pi |) is the average thermal momentum of the ath species of hadron. Ua is the number density of a*^ hadronic 
species and n is the total number desnity with n = Fig. |4Kb) shows shear viscosity to entropy density ratio 

estimated in two models, EHRG-I and EHRG-H. We note that the effect T (and p) dependent hadron masses is also 
reflected in transport properties. Shear viscosity is proportional to average thermal momentum which is certainly 
affected by temperature dependent hadronic species in the system. Nevertheless, the shear viscosity itself does not 
change much [Fig. UKa)], but the ratio rj/s is smaller in EHRG-H model than in EHRG-I due to more rapid increase in 
the entropy density in former [Fig. |2l[c)]. This effect is more important around transition temperature since the shear 
viscosity shows peculiar behavior around this temperature. It may be interesting to compare these results with the 
results of Ref.[l6l| where the authors estimated r]/s within EHRG model extended by inclusion of exponentially rising 
Hagedorn density of states. They observed that the inclusion of Hagedorn density of states significantly lowers 77 /s 


^ Alternating Gradient Synchrotron 
® Super Proton Synchrotron 






7 


and this ratio approaches close to the KSS bound near Tc- Thus, they argued that the inclusion of Hagedorn states 
could explain the low value of shear viscosity in the hadronic phase. Since we observed the same behavior of 77/5 near 
Tc but with the inclusion of medium dependent hadronic states, it may be again tempting to conclude that the effects 
of Hagedorn states can be alternatively simulated by including T and ^ dependent hadron masses in EHRG. 




T(GeV) T(GeV) 

(a) (b) 


FIG. 4: (Color online)Left panel shows shear viscosity as a function of temperatnre in EHRG I and EHRG H. Right panel 
shows shear viscosity to entropy density in two models. 









IV. SUMMARY AND CONCLUSION 


In present work we studied the effect of temperature (T) and baryon chemical potential (p) dependent hadron 
masses on the thermodynamics of hadronic matter. We used SU(3) Nambu-Jona-Lasinio model to compute T (and 
/i) dependent masses of constituent quarks (u, d, s) as well as light mesons (tt, K, rj and rj' ). For heavier mesons 
and baryons we used linear scaling rule in terms of constituent quarks. We found that although the general behavior 
don’t change much, taking into account T (and fi) dependent hadron masses makes the thermodynamical quantities 
numerically larger than with zero temperature (and chemical potential) hadron masses at moderately high temper¬ 
atures. This effect is more pronounced above T= 0.14 GeV. We argued that this behavior can be explained by 
considering the Boltzmann factor Exp(-M/j/T) which is the measure of probability of thermal excitation of specific 
hadronic species at given temperature. Further, we found that including the medium effects of hadrons in EHRG 
explains LQCD data qnite well over wide ranges of temperatures. Since it has been observed that the supression 
effect in EHRG together with the enhancement effect due to Hagedorn states leads to better agreement with LQCD, 
it may be tempting to conclude that the effects of Hagedorn states can be alternatively simulated by including T and 
dependent hadron masses in EHRG. At finite chemical potential, the EHRG-H is in better agreement with LQCD 
simulations than EHRG-I. We further emphasized the importance of taking into account this effect while estimating 
transport properties of hadronic matter especially around transition temperature. 
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